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Abstract 

The T2K neutrino oscillation experiment will start in 2009. In this experiment 
the center of the neutrino beam from J-PARC at Tokai village will go through 
underground beneath Super-Kamiokande, reach the sea level east of Korean shore, 
and an off-axis beam at 0.5° to 1.0° can be observed in Korea. We study physics 
impacts of putting a 100 kt-level Water Cerenkov detector in Korea during the 
T2K experimental period. For a combination of the 3° off-axis beam at SK with 
baseline length L = 295km and the 0.5° off-axis beam in the east coast of Korea 
at L = 1000km, we find that the neutrino mass hierarchy (the sign of — mf ) 
can be resolved and the CP phase of the MNS unitary matrix can be constrained 
uniquely at 3-0" level when sin^ 26'rct ~ 0.06. 
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The results of solar and atmospheric neutrino oscillation experiments are consistent 
with the 3 neutrino model, so are all the other observations except for the LSND ex- 
periments pp. In this paper we assume the 3 neutrino model, which has 6 observable 
parameters in neutrino oscillation experiments. They are 2 mass squared differences, 
3 mixing angles and one CP phase. The atmospheric neutrino oscillation experiments 
which measure the i/^ survival probability determine the absolute value of one of the two 
mass squared- differences and one mixing angle P| as 

1.5 X 10"^ <\ml-ml\< 3.4 x 10"^ and sin^ 2^atm > 0.92 (1) 

at the 90% confidence level. The K2K experiment |3J, which is the first accelerator 
based long baseline (LBL) neutrino oscillation experiment, confirms the above results. 
The solar neutrino experiments, which measure the z/g survival probability in the sun 4], 
and the KamLAND experiment which measures the Ue survival probability from reactors 
determine the other mass squared- difference and another mixing angle as 

ml-ml = 8.2+° j x IQ-^eV^ and tan^ = 0.401°;°? . (2) 

The CHOOZ reactor experiment gives the upper bound of the third mixing angle 
(6'rct) as 

sin^ 2^rct < 0.20 for - = 2.0 x 10"^eV^ , 

sin^ 20,ct < 0.16 for \ml - = 2.5 x lO'^eV^ , (3) 

sin^ 29,^t < 0.14 for \ml - m?| = 3.0 x 10" , 

at the 90% confidence level. The CP phase (^mns) been constrained. In the 

future neutrino oscillation experiments, we should not only measure sin^ 26*^1 and S^^^ , 
but also resolve the parameter degeneracies [111313, ^^ch as the sign of m| — mf. 

There are many next generation LBL neutrino oscillation experiments ^U]; which 
plan to measure the model parameters and to solve the parameter degeneracy. In this 
paper, we investigate the possibility of detecting in Korea the neutrino beam from J- 
PARC (Japan Proton Accelerator Research Complex) at Tokai village ^1], that will 
be available during the period of the T2K (Tokai-to-Kamioka) experiment ^21- In the 
T2K experiment, the center of the neutrino beam will go through underground beneath 
Super-Kamiokande, and reach the sea level near the Korean shore. At the baseline length 
L = 295km away from J-PARC, the upper side of the beam at 2° to 3° off-axis angle 
is observed at Super-Kamiokande (SK), and the lower side of the same beam at 0.5° 
to 1.0° appears in the east coast of Korea ^H]; at about L = 1000km; see Fig. 1. In 
order to quantify our investigation, we study physics impacts of putting a lOOkt water 
Cerenkov detector in Korea during the T2K experiment period, which is for 5 years with 
10^^ POT (protons on target) per year. 

The merits of measuring the T2K beam in Korea can be summarized as follows: 

1. Because 0.5° — 1.0° off-axis beam has significantly harder energy spectrum than 
2.5° — 3.0° off-axis beam, one can measure the i/^ transition probability near 
the oscillation maximum both at Korea and at SK, at the same time. 

2. Because of the matter effect that grows with the baseline length, the difference in 
the i'fj_ — > Ue transition probability between Korea and SK can reveal the neutrino 
mass hierarchy pattern [T3] . 
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Figure 1: The off-axis angle of the neutrino beam from J- PARC on the sea level in 
Korea, when the beam center is 2.5° (left) and 3.0° (right) off at SK. The baseline length 
for L = 1000, 1100, 1200km are shown by vertical contours, and the off-axis angles are 
shown by elliptic contours between 0.5° and 4.0°. 

3. The energy dependence of the oscillation probabilities measured by selecting the 
quasi-elastic events both at Korea and at SK allows us to constrain both cosine 
and sine of the CP phase ^^ns- 

The z/^ survival probability and the z/^ Ue transition rate can be expressed as 

P.,^u, = 1 - sin2 2^,tm (1 + A^) sin2 + B>^^ , (4) 

P. = 4 sin^ e^tn. sin^ 6,,^ {I + A') sin^ + bA , (5) 



where Aj^ = (m^ — mf)L/2E. Here A"" and B°' are the corrections to the magnitude 
and the phase of the oscillation probabilities, respectively, from the matter effect and the 
smaller mass-squared difference. If we keep only those terms linear in A12 and aL/2E, 
we find 

aL 1 - 2sin2 6'atm . 2 a 
A^^ = sm fe'rct , 

COS^ 6'atm 

< B^ = ^ l-2sin ^atm ^.^2 g^^^ (6) 

4^ C0S2 6'atm 

^ (cos^ 6*801 + tan^ 6'atm sin^ 6'soi sin^ 6'i.ct - tan 6'atm sin 26'soi sin 6'rct cos 6^ 



"mns 



Ae n o • 2/) N ^12 sin2^sol . . 

Ai3£/ Z tan t^atm sm t^rct 



I ^iii-^ ^ L.an (7atm sm f/j-ct /y\ 

. ^ =-4E^^-2^^^ ^-^^ + — (, 2tangatmsing..t '"^^-""" 



^MNs ~ sin 6'soi 
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The angles are expressed as the terms of the 3x3 MNS matrix elements sin^ O^ct = 
|t/e3p, sin^ 6'atm = I^^msP, sin^ 26'soi = 4|[/ei^^e2p as in refs.[T71[THl. A13 > for the normal 
hierarchy, A13 < for the inverted hierarchy, and A12 ~ |Ai3|/30 from the constraints 
eqs. (jH) and 0. The term a controles the matter effect [T^ . 



a = 2V2GFEne = 7.56 x IQ-^eV^ 



^g/cmO (cev) 



(8) 



where is the number density of the electron and p is the matter density. The magnitude 
of aL/2E is about 0.17 at SK, while it is about 0.57 at Korea with p = 3.0g/cm'^. By 
inserting the typical values of the observed parameters eqs (0) and Q, we find 
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near the oscillation maximum, IA13I ~ vr. Because the magnitude of A^^ and are 
very small, the z/^ survival rate is rather insensitive to the matter effect and subleading 
terms of order A12. Accordingly, measurements of the survival probability improve 
the constraints on |m| — and sin^ 26'atm- On the other hand, both A^ and B^ can 
affect the — > z/g oscillation probability significantly. Most importantly, the magnitude 
of the transition probability receives the correction term from the matter effect whose 
sign follows that of mg — ml and whose magnitude grows with L near the oscillation 
maximum, IA13I ~ tt. If we define the difference of the v^. transition probability 

between at Korea and SK, it can be expressed near the oscillation maximum as 



APnormal — -P/i^e(-^^far, A13 — +7r) — P^^e(Lncar7 A13 — +7r) , 
APinvcrted = P ii.-^eij-'iw^ A13 = — Tt) — P^^e(-^ncar, A13 = — Tt) , 



respectively, for the normal hierarchy (A13 = tt) and the inverted one (A13 = — vr). The 
difference is positive, and can be expressed as 



APnormal ~ APinvertcd ~ 8 siu^ 6'atm SIH^ ^rct 




(^w^ ) ^ "'o^n.r^^ ) • (12) 
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The difference grows linearly with the distance, Lfar, as long as the oscillation maximum 
is covered by the ffux. The ability of excluding the inverted hierarchy (A13 = — tt) is 
then determined by the error of the APnormah which can be estimated as 



1/2 



P^i^e (-^near ) 



1/2 



(13) 



Here N^. is the number of z/g appearance events. N]^^ /N'^'^^'^ can be expressed as 
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(14) 



where V denotes the fiducial volume of the detector and ^{E^, L) is the neutrino beam 
flux at L, which is proportional to (1/L)^. The neutrino cross section of Quasi-Elastic 
events is almost the constant in the 0.7 - 10 GeV region. Typical event number, N^^'^'^, 
for sin^ 2^rct = 0.1 and 6^^^^ = 0° during the 5 years running at SK is about 200; see a 
few bins around 0.5 GeV in Fig. (j2I). We therefore estimate significance of excluding the 
fake hierarchy as 
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(15) 

We find that when we put a 100 kt detector at L = 1000 km, the significance can exceed 
3.5 in this rough estimate. 

As of March 2005, there is no proposal to construct a huge neutrino detector in 
Korea. In our case study, we consider a 100 kt level detector, in order to compensate for 
the decrease in the neutrino flux which is about (300 km/1, 000 km)^ ~ 1/10 of that at 
SK. We adopt a Water Cerenkov detector because it allows us to distinguish clearly the 

events from events. We use the Charged-Current-Quasi-Elastic (CCQE) events in 
our analysis, because they allow us to reconstruct the neutrino energy by measuring the 
strength and the orientation of the Cerenkov lights. Since the Fermi-motion effect of the 
target dominates the uncertainty of the neutrino energy reconstruction, which is about 
80 MeV, in the following analysis we take the width of the energy bin as 5Ey = 200 MeV 
for > 400 MeV. The signals in the i-th energy bin. El = 200MeV x (i + 1) < E^ < 
El + 6 El,, are then calculated as 



Nl{v,) = MNa j^^ ^.,{E) P.^^.SE) a^'^iE) dE , 



(16) 



where Py^^^^ is the neutrino oscillation probability including the matter effect, M is 
the mass of detector, A^i = 6.017 x lO^'^ is Avogadro constant, is the t'^ flux from 
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Figure 2: The typical numbers of the /i events (left), and those of the e events (right), 
for the exposure time of 5 years, for the 3.0° OAB at SK (open square), and for the 0.5° 



OAB at L 
and = 



1000 km with lOOkt detector (solid circles). The inputs are sin 26'i.ct = 0.1 



0°. 



J-PARC, and is the CCQE cross section per nucleon in water. In this study, we 
assume that the fiducial volume of Super-Kamiokande is 22.5 kt, and that of a detector 
in Korea is 100 kt. We also assume that the detection efficiencies of both detectors for 
the CCQE events is 100%. The typical event numbers with sin^ 26'rct = 0.1 and (5^j^g = 0° 
is shown by Fig. |21 

We consider the following background events for the signal e- and yU-like events 

iV^'^^ = Kiu,) + Nl{v,) + Nl{v,) , (a = e , /i) . (17) 

The three terms represent the contribution from the secondary neutrino flux of the 
primary beam, which are calculated as in eq. (fTB|) where $y^(i?) is replaced by (^^^{E) 
for Up = Ue ,Ue ,i^fi- After summing up these background events, the e-like and /i-like 
events for the i-th bin are obtained as 

K = Ki^,) + N'^'''' ^ (a = e,/i). (18) 

Since our concern is the possibility to distinguish the neutrino mass hierarchy and to 
measure sin^ 29rct and the CP phase, we study how the above 'data' can constrain the 
model parameters by using the function 

^X^ = xIk + XKv + Xsys + Xjara • (19) 

Here the first two terms, Xsk Xxr; measure the parameter dependence of the fit to 
the SK and the Korean detector data. 
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where the summation is over all bins from 0.4 GeV up to 5.0 GeV for A^"^, 1.2 GeV for 
Nf. at SK, and 2.8GeV for at Korea. These upper bounds are chosen such that most 
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of the bins used in our analysis contain more than 10 events. Here A'^ ^ is the calculated 
number of events in the i-th bin, and its square root gives the statistical error. In our 
analysis, we calculate Ai'^ g by assuming the following input ('true') parameters: 
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(21) 



with the constant matter density, p^'^'^^ = 2.8 g/cm^ for T2K and p^™*^ = 3.0 g/cm^ for 
the Tokai-to-Korea experiments. Note that in eq. (j2T| . we assume the normal hierarchy 
(777.3 ~ "^1 > 0) examine several input values of sin^ 26'rct and 6^,^^ . 

A^f* is calculated by allowing the model parameters to vary freely and by allowing 
for systematic errors. In our analysis, we consider 4 types of systematic errors. The 
first ones are for the overall normalization of each neutrino flux, for which we assign 3% 
errors, 

= 1±0.03, (22) 

for i>j3 = Ue, Ue, z/^, z/^, which are taken common for T2K and the Tokai-to-Korea ex- 
periments. The second systematic error is for the uncertainty in the matter density, for 
which we allow 3% overall uncertainty along the baseline, independently for T2K (/p^) 
and the Tokai-to-Korea experiment (f^^)'- 



Pf = fip'r (^ = SK, Kr). (23) 
The third uncertainty is for the CCQE cross section, 

^QE, fit ^^QE^QE, true _ (24) 

Since Ve and QE cross sections are expected to be very similar theoretically, we assign 
a common overall error of 3% for z/e and {ff^ = /^^ = ff^), and an independent 3% 

error for z/g and z/^ QE cross sections (/^^ = /^^ = /^^). The last one is the uncertainty 
of the fiducial volume, for which we assign 3% error independently for T2K (/y^) and 
the Tokai-to-Korea experiment (/y"^). A'"^''^* is then calculated as 



H^/3) = f^"" /f '^^ N^M , (25) 



and Xsys has four terms; 



Asys 



E 





+ 



7^ 



V 



0.03 



(26) 



7 



20 r 



16 



4 - 



OAB 2.5 at SK 

000km • 

050km o 

1QQkm A 

150km □ 



05 TD r5 215 Z5 315 

off-axis angle in Korea (degree) 



20 r 



16 



OAB 3.0 at SK 



000km 
050km 
km 
km 



05 TD r5 215 23 JIG 

off-axis angle in Korea (degree) 



Figure 3: Minimum Ax^ as functions of the off-axis angle and the base-hne length from 
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mg — rrii 



J-PARC at Tokai, when the normal hierarchy 
in generating the events, and the inverted hierarchy {m\ 
fit. The left-hand figure is for the 2.5° off-axis beam at SK, and the right-hand one is 
for the 3.0° beam. 



2.5 X 10"^eV > 0) is assumed 
: — ml < 0) is assumed in the 



In short, we assign 3% errors for the normalization of each neutrino flux, the z/g and 
Ue CCQE cross sections, the effective matter density along the base line, and for the 
fiducial volume of SK and the Korean detector. Finally, Xpara accounts for the present 
constraints on the model parameters, summarized in eqs. (1) and (2): 



A- 



|(m| — mf)^*| — \ {ml — m^)*™'^|\ / {ml — mf)*^* — {m^ — mlY^^^^ 



P^''' \ 0.5 X 10-3 \ 0.6 X 10-5 



^ / sin^ 2etL - sin^ 2^ \ ^ / sin^ 2^ - sin^ 2g \ ^^7) 



Here we interpret the 90% CL lower bound on sin^ 2^atm in eq. as the 1.96cr constraint 
from sin^ 2^atm is greater than 0.96, and the asymmetric error for tan^ ^soi in eq. (j2)) has 
been made more symmetric for sin^ 26'soi • We do not include the bounds on sin 2,6]^ct in 
eq. (jni) in our Ax^ function. In total, our Ax^ function depends on 16 parameters, the 
6 model parameters and the 10 normalization factors. 

First, we search for the best place for the detector in Korea and the best combination 
of the off- axis angle for SK and the Korean detector to determine the sign of mg — m^. 
We show in FiglHlthe minimum Ax^ as functions of the off-axis angle and the base-line 
length in Korea, when the data, A^^, are generated for the normal hierarchy, mg — mf = 
2.5 X lO'^eV^ > 0, eq. ()21|). and the fit has been performed by assuming the inverted 
hierarchy, m| — mf < 0. We set sin^2^J:™^ = 0.10 and = 0° in this analysis. The 
left hand figure shows the minimum Ax^ for the 2.5° off-axis beam at SK, and the right 
hand one is for the 3.0° off-axis beam at SK. The four symbols, solid circle, open circle, 
triangle, and square are for L = 1000km, 1050km, 1100km, and 1150km, respectively. 
When the off-axis angle at SK is 2.5°, the 0.5° beam does not reach the Korean coast; see 
Fig. ^ It is clear from these figures that the best discriminating power is obtained for 
the combination L = 1000km and 0.5°, which is available only when the off-axis angle at 
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Figure 4: Minimum Ax^ as a function of sin^ 26'^*^ when the normal hierarchy {m\ — m\ = 
2.5 X 10~^ eV'^ > 0) is assumed in generating the events, and the inverted hierarchy 
(m| — ml < 0) is assumed in the fit. The 4 hues are for 4 input CP phase values, 
5^rue ^ 0.0° (sohd), 90° (long-dashed), 180° (short-dashed) and 270° (dotted). The left 
figure is for sin^ '^^llt'^ = 0.10 and the right one is for sin^ 26'*^;^'^ = 0.06. 

SK is 3.0° (right figure). For this combination, we can distinguish the inverted hierarchy 
from the normal one at more than Aa level. These figures show that longer base-line 
gives larger minimum Ax^ for the same off- axis angle. This is because of the increase in 
the matter effect. For the same base-line length, lower off axis angle beams give better 
discriminating power. This is because the neutrino flux with smaller off-axis angle is 
harder fT^ 1^ . and the stronger matter effect to help us to distinguish the neutrino 
mass hierarchy fn\ ^ ■ 

Here after, we study the prospect for measuring the sign of 171"^ — mf in more detail for 
the best combination, L = 1000km and 0.5° in Korea, and 3.0° for SK. Fig. |3] shows the 
minimum Ax^ as functions of the fitting parameter sin^ '^^fct by assuming the inverted 
hierarchy, when the normal hierarchy (mg — ml = 2.5 x 10~^ eV'^ > 0) is assumed in 
calculating the event numbers. There are 4 lines in each figure, which correspond to 4 
input values of the CP phase, = 0° (solid), 90° (long-dashed), 180° (short-dashed) 
and 270° (dotted), respectively. The left figure is for sin^ 26'|:j;"'^ = 0.1 and the right one 
is for sin^ 2^*^"^ = 0.06. Ax^ is mainly controled by the difference of A^* between the 
normal hierarchy and the inverted hierarchy which is proportional to sin^ 6'rct; see eq. (I12|l . 
Because P^^^^^ in the inverted hierarchy is smaller than that in the normal hierarchy due 
to the matter effect, the fitting parameter, sin^2^fp\, tends to be larger than the input 
value, sin^26'*™'^. Since large sin^ 26'rct will be constrained more strongly in the future 
reactor experiments [^SIISI, we can conclude that the neutrino mass hierarchy will be 
determined at even higher confidence level once results from these reactor experiments 
are available. 

We also examine the capability of the Tokai-to-Korea LBL experiments for measuring 
the CP phase. We show in Fig. El regions allowed by this experiment in the plane of 
sin^26'i.ct and 5^i^g. The mean values of the inputs are calculated for the parameters of 
eq. (|2ip. In each figure, the input points (sin^2^J:™^, 5^^g) are shown by solid-circles for 
sin^26'J:™^ = 0.10, and 0.06. The regions where the minimum Ax^ is less than 1, 4, 9 are 
depicted by solid, dashed and dotted boundaries, respectively. Even though we allow the 
sign of m| — ml to vary in the fit, no solution with the inverted hierarchy that satisfy 



9 













'/' 


5 






( 












































110 






























































































80 
































































































90 
































































































o 



































0.02 


0.04 


0.06 


0.08 


0. 


10 


0.12 


0.14 



360 



270 



=1 80 



CO 



90 







0.02 0.04 0.06 o0-08 0.10 0.12 0.14 

sin^2act 



































































70 






























































































80 


























































































90 














































































































360 



270 



180 



90 



((OlKCO 



0.02 0.04 0.06 „0.08 0.10 0.12 0.14 
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at SK and Korea are calculated for the parameters of eq. (PT|). In each figure, the input 
points (sin^2^*™'^, 5^^^) are shown by solid-circles, and the regions where the minimum 
Ax^ is less than 1, 4, 9 are depicted by solid, dashed and dotted boundaries, respectively. 
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^Xmin < 9 appear in the figure. 

From these figures, we learn that S^^^^ can be constrained to ±30° at la level, when 
sin^ 26**™^ > 0.06. It is remarkable that we can constrain both sin^^^^g and cos(5^j^g with- 
out using anti-neutrino experiments. We can determine sin(5^j^g uniquely by measuring 
the z/g oscillation probability near the oscillation maximum both at SK and Korea. 
This is because the significant difference in the matter effect term in eq. (|7j) between SK 
and Korea allows us to resolve the correlation between sin^ 6'rct and sin5^j^g [211 123 12^1 • 
As for cos 5^f^g , it appears only in the phase shift of the z/^ z/g oscillation probability; 
see the term B'^ in eq. ((Tj). It is therefore important to measure the neutrino energy by 
CCQE events, in order to constrain cos 6^^^ . 

In this paper, we study the possibility of solving the degeneracy of the neutrino mass 
hierarchy and constraining sin^ 29,.^ and (5^j^g by measuring the T2K off-axis beam in 
Korea. We find that by placing a 100 kt level water Cerenkov detector in the east coast of 
Korea, we can determine the sign of ml — ml and constrain sin^ 26*1^ and S^j^^ uniquely, 
if sin^ 2^rct ^ 0.06 . 

Our results are based upon a very simple treatment of the systematic errors where 
3% overall errors are assigned for all the 10 normalization factors of eq. (j2SI)- Even if 
we enlarge all the systematic errors to 10% except for the matter density uncertainties, 
the significance of the mass hierarchy determination is not affected much. For instance, 
the Ax^in for the combination of the 3° OAB at SK and the 0.5° OAB at L = 1000 km 
in Fig. 121 is found to reduce from 18 to 16. Among the potentially serious background 
which we could not estimate in this paper are; 

• possible miss-identification of NC 7r° production as z/g CCQE events, 

• possible miss-identification of soft vr emission events as z/g CCQE events. 

Although the above uncertainties were found to be rather small at K2K experiments 
[27], we should expect them to be more serious at high energies. Dedicated studies of 
their effects on the neutrino-energy reconstruction efficiency are mandatory. In addition, 
careful studies including possible energy dependence of the flux and cross section uncer- 
tainties, location dependence of the matter density may be needed to justify the physics 
case of our proposal. 
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Note added: 
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